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One of the most fundamental phenomena and a reminder of the electron’s relativistic nature is 
the Rashba spin splitting for broken inversion symmetry. Usually this splitting is a tiny relativistic 
correction, hardly discernible in experiment. Interfacing a ferroelectric, BaTi 03 , and a heavy 5 d 
metal with a large spin-orbit coupling, Ba( 0 s,Ir) 03 , we show that giant Rashba spin splittings are 
indeed possible and even fully controllable by an external electric field. Based on density functional 
theory and a microscopic tight binding understanding, we conclude that the electric field is amplified 
and stored as a ferroelectric Ti-0 distortion which, through the network of oxygen octahedra, also 
induces a large Os-O distortion. The BaTiOs/BaOsCU heterostructure is hence the ideal test station 
for studying the fundamentals of the Rashba effect. It allows intriguing application such as the 
Datta-Das transistor to operate at room temperature. 

PACS numbers: 73.20.-r, 73.21.-b, 79.60.Jv 


An electric field control of the spin degree of freedom 
is the key to spintronics and magnetoelectrics [1]. In 
the prototypical spintronic device, the Datta-Das spin 
transistor [2], an electric field tunes the Rashba spin 
splitting and with that the spin precession frequency. 
The precession in turn controls the spin polarized cur¬ 
rent between two ferromagnetic leads. Microscopically, 
the Rashba spin splitting originates from the spin or¬ 
bit coupling (SOC) in a two dimensional electron gas 
(2DEG) with broken inversion symmetry perpendicu¬ 
lar to the 2DEG plane [3, 4]. It has been observed 
for metal surfaces [5, 6], semiconductor and oxide het¬ 
erostructures [7-12], and even in polar bulk materials 
[13]. The Rashba effect splits parabolic bands into two 
subbands with opposite spin and energy-momentum dis¬ 
persions E ± ( k) = (A 2 k 2 /2m*) ±<a#|k|. Here, m* is the 
effective mass, k the wave vector in the 2DEG plane, and 
aR the Rashba coefficient which depends on the strength 
of SOC and inversion asymmetry. An electric field mod¬ 
ulates this inversion asymmetry and consequently the 
Rashba spin splittings. The electric-field induced change 
of aR is however weak: up to 10 _2 eVA in semiconductor 
[7] or 3d transition metal oxide heterostructures [8-10]. A 
large, electric-field tunable Rashba effects are also much 
sought-after in the research area of topological insulators 
[14-16]. 

Giant Rashba effects with aR of the order of leVA have 
been reported for metal surfaces [6], Bi adlayers [17] and 
bulk polar materials BiTel [13, 18]. The large aR here 
relies on the surface or interface structural asymmetry, 
which is hardly changes in an external electric field. To 
enhance the tunability by an electric field, Di Sante et 
ai [19] hence suggested a ferroelectric semiconductor 
GeTe. For GeTe, an external electric field switches be¬ 
tween paraelectric and ferroelectric phase, breaking in¬ 


version symmetry and tuning on the Rashba spin split¬ 
ting. At first glance, this perfectly realizes an electric 
field control of a giant Rashba effect. However, there is 
an intrinsic difficulty: A single material cannot be both, 
a conductor with large Rashba spin splitting and a fer¬ 
roelectric which necessarily is insulating. 

In this letter, we propose to realize a giant switch- 
able Rashba effect by heterostructures sandwiching a thin 
metallic film of heavy elements in-between ferroelectric 
insulators, see Fig. 1. The thin film provides a 2DEG 
with strong SOC, while the inversion asymmetry is in¬ 
duced by the structural distortion of the ferroelectrics 
and tunable by an electric field. As a prime example, 
we study a heterostructure of transition metals oxides, 
BaOsOs/BaTiOs. BaTiOs is a well-established ferro¬ 
electric [20] . It has a simple high-temperature perovskite 
structure, with Ba atoms at the edges, a Ti atom at the 
center and the O atoms at the faces of a cube. Such 
a structure has an inversion symmetry center at the Ti 
site. At room temperature, inversion symmetry is bro¬ 
ken since a ferroelectric structural distortion occurs with 
a sizable Ti-0 displacement ZTi-o along one of the the 
cubic axes. BaOsOs has been recently synthesized; it is 
a metallic perovskite with four Os 5 d electrons and no 
sign of magnetism [21]. It has a perfect lattice match 
with BaTiOs. Both materials have the same cation, Ba, 
which will substantially reduce disorder during epitax¬ 
ial growth. For BaOsOs/BaTiOs we find an electric-field 
tunable Rashba spin splitting which is at least one magni¬ 
tude larger than the current experimental record [7-10]. 
The mechanism behind is nontrivial and summarized in 
Fig. 1. Substituting BaOsOs by BalrOs and BaRuOs 
yields a similar effect; the heterostructure can also be 
further engineered by varying its thickness and strain. 

Method. We mainly focus on a BaOsOs/BaTiOs het- 
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FIG. 1: Schematics of the mechanism behind the giant switch- 
able Rashba effect: (1) An external electric field is amplified 
and stored as a ferroelectric distortion, in particular Ti-0 dis¬ 
placements. (2) These displacements entail, via the oxygen 
octahedron network of the perovskite heterostructure, Os-O 
displacements. (3) As a consequence of the Os-O displace¬ 
ments the Os orbitals deform, which reflects the broken in¬ 
version symmetry. This orbital deformation and the strong 
SOC of Os finally lead to a giant Rashba spin splittings that 
is controllable by an electric-field. 


erostructure which has two structural phases: a paraelec- 
tric perovskite phase and a distorted ferroelectric phase, 
see Fig. 2(c) and (d), respectively. We use density- 
functional theory (DFT) with generalized gradient ap¬ 
proximation (GGA) potential [22] in the Vienna Ab ini¬ 
tio Simulation Package (VASP) [23, 24] and fully relax 
all the atomic positions, only fixing the paraelectric or 
ferroelectric symmetries. At zero temperature, the ferro¬ 
electric distorted perovskite has the lower GGA energy, 
but -as in the bulk- the undistorted paraelectric phase 
will prevail at elevated temperatures. 

Based on the fully relaxed atomic structures, electronic 
band structures are calculated with modified Becke- 
Johnson (mBJ) [25] exchange potential as implemented 
in the Wien2k code [26], which improves the calculated 
bandgap of BaTiOs. The SOC is included as a pertur¬ 
bation using the scalar-relativistic eigenfunctions of the 
valence states. Employing wien2wannier [27], we project 
the Wien2k bandstructure onto maximally localized [28] 
Wannier orbitals, from which the orbital deformation is 
analyzed and a realistic tight binding model is derived 
[29]. We also vary the thickness of the thin films, re¬ 
place the BaOsOs by BaRuOs or BalrOs, and simulate 
the strain effect by fixing the in plane lattice constant to 
the value of a SrTiOs substrate. Correlation effects from 
local Coulomb interaction are studied in the Supplemen¬ 
tary Material. 

Paraelectric phase. We mainly fo¬ 

cus on a BaOsOs/BaTiOs heterostructure, 
(Ba 0 s 03 )i/(BaTi 03 ) 4 , which consists of one BaOsOs 
layer alternating with four BaTiOs layers. Bulk BaTiOs 


is a dP insulator: the empty Ti 3d states lie about 3eV 
above occupied O 2 p states. Bulk BaOsOs is a d A metal: 
four electrons in the Os 5 d orbitals of t^g character are 
0.8eV above the filled O 2 p states. Since BaOsOs and 
BaTiOs share oxygen atoms at the interface, the O 2 p 
states align; and the Os 5 d states will stay in the energy 
gap of BaTiOs. Hence, the density functional theory 
(DFT) calculated band structures of BaOsOs/BaTiOs 
in Fig. 2 shows three Os t^g (■ xy , yz, xz) bands near the 
Fermi level; they are dispersionless along the z direction 
(not shown). That is, the low energy electronic degrees 
of freedom are confined by the insulating BaTiOs layers 
to the 0s0 2 plane, forming a 2DEG. Already in the high 
temperature paraelectric phase, this heterostructure 
confinement reduces the initial cubic symmetry of the 
Os t2 g orbitals in the bulk perovskite: At V the xy band 
in Fig. 2 (a) is 1.2eV lower in energy than the degenerate 
yz/xz bands (Fig. 2(a)). Including the SOC, the yz/xz 
doublet splits into two subbands with mixed yz/xz 
orbital character, see Fig. 2(b). For the paraelectric 
phase of Fig. 2(a-c), the spin degeneracy is however still 
preserved. 

To better understand the DFT results, we construct 
an interface hopping Hamiltonian Hq [30] for the 2DEG 
confined Os t^g electrons by a Wannier projection. 
The energy-momentum dispersion for the xy orbital is 
e(k ) xy = — 2 ti cos k x — 2t\ cos k y — Ats cos k x cos k y , while 
that of the yz is e(k) yz = — 2 t 2 cos k x — 2t\ cos k y (e(k) xz 
is symmetrically related by x y). The largest hop¬ 
ping t\ = 0.392eV is along the direction(s) of the or¬ 
bital lobes, t 2 = 0.033eV and £3 = 0.094eV indicate a 
much smaller hopping perpendicular to the lobes and 
along (1,1,0), respectively. The xy-yz/xz orbital split¬ 
ting at T is 2ti — 21 2 + 4^3 = l.leV and arises from the 
anisotropy of the orbitals. Overall, the three parameter 
tight binding model in Fig. 3(a) (dashed line) is consis¬ 
tent with DFT in Fig. 2(a). Next, we include the atomic 
SOC Hamiltonian • s, expressed in the t^g basis 

with £ = 0.44eV for Os. It lifts the yz/xz degeneracy 
and yields good agreement with DFT, cf. Fig. 3(a) (solid 
line) and Fig. 2(b). 

Ferroelectric phase. In sharp contrast to the para¬ 
electric case, ferroelectrically distorted BaOsOs/BaTiOs 
exhibits evident spin splittings in the presence of SOC, 
see Fig. 2 (f). Note, in the absence of SOC the DFT 
band structure in Fig. 2(e) is still very similar to the 
paraelectric case in Fig. 2(a). The spin splitting E^ 
of the xy band is small around T, but strongly en¬ 
hanced up to 50meV around the xy/yz crossing region 
as shown in Fig. 2(g), cf. Table I, . The other two sub¬ 
bands of yz/xz character also show a strong splitting 
around T. In this respect, the upper yz/xz subband 
exhibits a standard Rashba behavior, see the magnifi¬ 
cation in Fig. 2(h), where the momentum offset fco = 
0.043 A -1 , the Rashba energy Fr = 4meV, and hence 
<aR = 2 Fr/&o = 0.186 eVA. The behavior of the lower 


















FIG. 2: DFT calculated band structures of (Ba 0 s 03 )i/(BaTi 03)4 multilayers for the paraelectric (c) and ferroelectric structure 
(d). The paraelectric heterostructure breaks cubic but not inversion symmetry so that the Os xy orbital splits off from two 
degenerate yz/xz orbitals without SOC (a); including SOC further lifts the yz/xz degeneracy (b). The ferroelectric state 
further breaks inversion symmetry, which essentially does not change the bandstructure in the absence of SOC (e) but leads to 
substantial spin splittings in the presence of SOC (f). The spin splitting E ^ at the xy — yz crossing region is magnified in (g), 
while the standard Rashba spin splitting Er with momentum offset ko around T of the upper yz/xz mixed orbital is magnified 
m fhV 




FIG. 3: Tight binding bandstructure of the three Os t 2 g or¬ 
bitals for (a) the paraelectric case with and without SOC and 
(b) the ferroelectric case with SOC and asymmetry parameter 
7 = 0.016 eV taken from the Wannier projection. 


yz/xz subband is more complex (distorted) due to the 
proximity of the crossing region. 

These spin splittings arise from ferroelectric distortions 
which break the inversion symmetry. Bulk BaOsOs has 
no such ferroelectric distortions and is inversion sym¬ 
metric with Os and O atom in the same plane. Also 
for the paraelectric heterostructure, the 0 s 02 layer is 
still an inversion plane. The ferroelectric heterostructure 
breaks inversion symmetry. As listed in Table I, the av¬ 
eraged ferroelectric Ti-0 displacement in BaTiOs layers 
is around 0.14 A, which efficiently induces a sizable Os-O 
displacement of 0.05 A via the oxygen octahedron net¬ 
work of the perovskite structure. This structural distor¬ 
tion modifies the crystal environment of Os t 2g orbitals 
and deforms the orbital lobes, see Fig. 1. 

In order to quantify the orbital deformation, we project 
the DFT results above onto maximally localized Wan¬ 
nier orbitals and directly extract the key parameter: a 
directional, spin-independent inter-orbital hopping 7 = 
(xy\H\yz(R)), were R is the nearest neighbor in x direc¬ 


tion. In k-space one gets as a matrix for the yz,xz,xy 
orbitals (independent of spin) [29] 

( 0 0 2 i sin k x 

0 0 2 i sin k y 

—2 ismk x —2ismk y 0 

We obtain 7 = 0.016 eV for the ferroelectric het¬ 
erostructure, while for the paraelectric case 7 = 0 due 
to inversion symmetry. Let us emphasize that 7 is the 
relevant measure for the orbital deformation and inver¬ 
sion symmetry breaking. In combination with the strong 
SOC of OS 5 d electrons, H 1 results in the present giant 
Rashba spin splitting. 

The Hamiltonian Hq + + H 1 already well describes 

the spin splittings of the xy orbital in good agreement 
with DFT results, but underestimates the spin splittings 
of the yz/xz subbands near T. The reason for this is that 
the atomic SOC is not an accurate description anymore 
because the SOC of Os is too strong and the deviation 
from the inversion symmetry too large. Hence, we need 
to go beyond the purely atomic SOC and include an in¬ 
version asymmetry correction to the SOC matrix so that 
H $ in the t 2g basis (yz | t), yz | |), xz\ t), xz\ |), xy | t), 
xy\\)) reads 



( 0 

2 y' sin k y 

-ii 

b 2 

0 
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£ 
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27 ' sin k y 
0 
0 

2 
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0 

0 

— 27 ' sin k x 



0 0 \ 

2 0 

—27'sinfca; 0 if 

0 if 0 

-if 0 0 

0 0 0 / 


For the bandstructure of Fig. 3(b) we have taken 7 ' = 
0.022 eV from the Wannier projection, yielding good 
agreement with the DFT results Fig. 2(f). 






















































4 


TABLE I: Polar distortions and Rashba spin splittings of (Ba 0 s 03 ) n /(BaTi 03 ) m multilayers, alternating n layers of BaOsOs 
and m layers of BaTiOs. Second and third column: averaged displacement of Os-O and Ti-0 along z as calculated by DFT. 
Forth column: strength of the interface asymmetry term 7 obtained from the Wannier projection. Fifth to seventh column: 
DFT calculated parameters characterizing Rashba effects: momentum offset ko, Rashba energy Fr, and Rashba coefficient cm 
of the upper yz/xz subband around T as defined in Fig. 2(h). Eighth column: spin splitting around the crossing region 
of xy and yz orbitals as in Fig. 1(g). To demonstrate the effects of heterostructure engineering, we also list the results of 
BaOsOs/BaTiOs strained by a SrTiC >3 substrate, as well as BaRuOs/BaTiOs and BalrOs/BaTiOs. Band structures are given 
in Supplementary Material. 


n : m 

Os-O (A) 

Ti -0 (A) 

7 (eV) 

fco (A- 1 ) 

E n (eV) 

<a R (eVA) 

El v (eV) 

1:3 BaOsOs/BaTiOs 

0.038 

0.065 

0.011 

0.025 

0.001 

0.08 

0.031 

1:4 

0.049 

0.103 

0.016 

0.043 

0.004 

0.186 

0.053 

2:3 

< 0.010 

< 0.010 

0.003 

0 

0 

- 

0.004 

2:4 

0.024 

0.071 

0.010 

0.035 

0.002 

0.114 

0.018 

1:3 strained 

0.095 

0.223 

0.030 

0.105 

0.021 

0.396 

0.102 

1:4 strained 

0.099 

0.229 

0.031 

0.115 

0.024 

0.417 

0.106 

1:4 BalrOs/BaTiOs 

0.115 

0.119 

0.021 

0.145 

0.053 

0.731 

0.051 

1:4 BaRuOs/BaTiOs 

0.082 

0.119 

0.015 

0.128 

0.016 

0.250 

0.007 


Electric-field tunability. The Rashba spin splittings 
can be tuned by an electric field because of the ferroelec- 
tricity of BaTiOs. It is a nature of ferroelectrics that an 
external electric field can tune and even reverse the polar¬ 
ization and structural distortion [17, 20, 31-33], including 
that of the BaOsOs layer. Because of the BaOsOs layer, 
the ferroelectricity of the heterostructure is not as strong 
as in bulk BaTiOs. Hence a smaller electric field change 
than in the bulk (E=10 3 V/cm [33]) is be needed for going 
through the ferroelectric hysteresis loop. Such an elec¬ 
tric field can tune the Rashba coefficient listed in Table 
I, which are at least one order of magnitude more larger 
than in semiconductors [7] or 3d oxide heterostructures 
[8-10]. For applications it is of particular importance 
that Er now clearly exceeds 0.025 eV, the thermal energy 
at room temperature. Furthermore, we propose to take 
advantage of the ferroelectric hysteresis and the remnant 
ferroelectric polarization, using our heterostructure as a 
spintronics memory device. 

Heterostructure engineering. Besides applying an elec¬ 
tric field, heterostructure engineering such as varying its 
thickness, strain, and the material combination is a pow¬ 
erful way to tune the Rashba spin splittings. First, we 
vary the thickness of BaOsOs and BaTiOs. It is well 
known experimentally and theoretically [20, 34] that be¬ 
low a few nanometer thickness of the BaTiOs film, the Ti- 
O displacement is reduced and eventually ferroelectric¬ 
ity vanishes. As shown in Table I, when decreasing the 
thickness of BaTiOs and increasing that of BaOsOs, the 
Rashba spin splittings is indeed substantially reduced. 
Second, ferroelectricity of BaTiOs is sensitive to strain 
[35]. Taking SrTiOs as a substrate provides a 2.5% in¬ 
plane compressive strain. This enhances the ferroelectric 
distortion as well as Rashba spin splitting, am=0.417eVA 
in Table I. Given the giant Rashba splitting, photoemis¬ 
sion and transport measurements should be able to val¬ 
idate this effect by comparing thin films with different 


thickness and substrate. Third, we replace Os by other 
5 d elements for shifting the Fermi energy; and also since 
some Os oxides are toxic. Indeed, similar physics can 
be found in BalrOs/BaTiOs with cm = 0.731eV A, see 
Table I, and in BaRuOs/BaTiOs. The latter also of¬ 
fers a platform to study the interplay of Rashba physics 
and electronic correlations; for instance, novel phenom¬ 
ena is expected in mixed 3d-bd transition-metal materials 
[36]. Note a strong temperature dependence of the qlr in 
SrlrOs was recently reported [37]. 

Conclusion. We propose to combine a heavy 5 d 
compound such as BaOsOs and a ferroelectric such as 
BaTiOs in a heterostructure for realizing giant Rashba 
spin splittings that are switchable by an electric field. 
The physics behind is that ferroelectric BaTiOs ampli¬ 
fies the electric field through a distortion (polarization) 
of its Ti06 octahedra. This distortion very efficiently 
propagates to the Os06 octahedron via the shared oxy¬ 
gens, breaking inversion symmetry and deforming the Os 
orbital lobes. Together with the strong SOC it leads to a 
giant tunable Rashba effect with, e.g., cm=0.731eVA for 
BalrOs/BaTiOs. This is at least one magnitude larger 
than state-of-the-art. Moreover, the ferroelectric hystere¬ 
sis and remnant polarization leads to a memory effect 
which can be exploited for spintronics memory devices. 
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